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Results of a numerical study are presented for hypersonic low-density flow about a 70-deg blunt cone using 
direct simulation Monte Carlo (DSMC) and Navier-Stokes calculations. Particular emphasis is given to the effects 
of chemistry on the near-wake structure and on the surface quantities and the comparison of the DSMC results 
with the Navier-Stokes calculations. The flow conditions simulated are those experienced by a space vehicle at an 
altitude of 85 km and a velocity of 7 km/s during Earth entry. A steady vortex forms in the near wake for these 
freestream conditions for both chemically reactive and nonreactive air gas models. The size (axial length) of the 
vortex for the reactive air calculations is 25% larger than that of the nonreactive air calculations. The forebody 
surface quantities are less sensitive to the chemistry than the base surface quantities. The presence of the afterbody 
has no effect on the forebody flow structure or the surface quantities. The comparisons of DSMC and Navier-Stokes 
calculations show good agreement for the wake structure and the forebody surface quantities. 


Nomenclature 

A = base area of cone, nd 2 /4 

C D - drag coefficient, 2D/p 00 V^ 0 A 

Cf = skin-friction coefficient, 2x w /p 00 

Ch = heat transfer coefficient, 2 q/pooV^ 

C p = pressure coefficient, 2p/p 00 V^ 0 

D = drag 

d = base diameter, 2 m 

d ie { = molecular diameter at reference temperature 

Kn = Knudsen number, X/d 

M = Mach number 

M = molecular weight of air, 28.96 g/mole 

N = Avogadro’s number, 6.02252 x 10 26 particles/kg-mole 

N = atomic nitrogen 

N 2 = molecular nitrogen 

O = atomic oxygen 

0 2 = molecular oxygen 

p = pressure 

q - heat flux 

Rb = cone base radius 

R c = corner radius 

Re = Reynolds number, pVd/p. 

Re 2 ~ total Reynolds number, Re 0 o(Poo/Pstag) 

R n = cone nose radius 

1Z = universal gas co nstant, 8.3 143 J/mole-K 

S = speed ratio, V y / M /21ZT 

s = distance along the body surface measured from the 
stagnation point 

s = temperature exponent of the coefficient of viscosity 

T - thermodynamic temperature 
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T ov = overall kinetic temperature 

T w = surface temperature 

u = axial velocity 

V = velocity 

v = radial velocity 

X = mole fraction 

x = axial distance from stagnation point measured along 
symmetry axis 

x a = location of rear stagnation point 

x b = location of wake stagnation point 

y = radial distance from symmetry axis 

T = Gamma function 

A = size of the vortex, (x b — x a )/d 

X = mean free path 

p = dynamic viscosity 

p = density 

cr = collision cross section 

r = shear stress 

Subscripts 

ref = reference value 

stag = stagnation point 

w = surface values 

oo = freestream values 

Introduction 

T HERE is a lack of experimental and computational data for 
blunt-body wake flows under rarefied flow conditions. Precise 
determination of wake closure is a critical issue for aerobrakes be- 
cause the low lift-to-drag ratio aeroshell designs impose constraints 
on payload configuration. The payload must fit into the wake cone 
to minimize the heating since a heating spike is generally asso- 
ciated with reattachment of the separated flow in the near wake. 
Because of the complicated nature of the wake flow a perception 
exists 1 that the wake aerothermodynamics cannot be predicted ac- 
curately. A number of fundamental issues exist concerning such 
flows: 1) What role does thermochemical nonequilibrium play in 
the near wake? 2) How does the wake structure change as a func- 
tion of rarefaction? 3) To what limits is the continuum modeling 
valid as rarefaction in the wake is progressively increased? Only re- 
cently, computational investigations 2-11 of blunt-body wake struc- 
ture have begun to answer some of the previous questions and to 
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isolate critical features of the wake structure of planetary probes and 
aerobrakes. 

Generally two methods are used to simulate the wake flows at high 
altitudes. One is the continuum approach in which a set of model 
equations is solved numerically, and normally this set is the Navier- 
Stokes equations. The Navier-Stokes modeling becomes inadequate 
for large local Knudsen numbers such as those occurring when the 
relatively high-density forebody flow expands into the wake. The 
second method is of a molecular nature where a direct physical sim- 
ulation is done by following the motion and interaction of modeled 
molecules. This method is equally valid at low and high local Knud- 
sen numbers, and therefore it is more appropriate to simulate the 
wake flows with the molecular method under rarefied conditions. In 
the present study, computational codes representative of both meth- 
ods are used to simulate hypersonic wake flows about a 70-deg blunt 
cone. The molecular simulations are achieved with Bird’s 12-14 direct 
simulation Monte Carlo (DSMC) method, whereas the continuum 
simulations are performed using the Navier-Stokes code developed 
by Olynick and Hassan 15 and Olynick et al, 16 

The purpose of the present paper is to provide an improved under- 
standing of the effects of chemistry on blunt-body wake structure 
at flight conditions and to compare the results with Navier-Stokes 
calculations with slip boundary conditions. In addition, the DSMC 
calculations are also performed for a 70-deg blunt cone with a cylin- 
drical as well as a proposed Mars Environmental Survey (MESUR) 
Pathfinder probe afterbody configurations. 

Computational Method and Boundary Conditions 

Both the molecular and continuum methods used in the present 
calculations are briefly described next. 

Molecular Method 

The direct simulation Monte Carlo (DSMC) method 12-14 is used 
for the present calculations. The method and requirements for appli- 
cation of DSMC have been presented in previous publications 17,18 
and are not repeated here. 

The computational domain used for the calculations is large 
enough to capture most of the body disturbance at the upstream 
and side boundaries. Thus, freestream conditions are specified at 
these boundaries. The flow at the downstream outflow boundary is 
supersonic, and vacuum conditions are specified. 

The flowfield was divided into 1 8 regions and a fine grid resolution 
was used for cells in the wake regions. The cell dimension normal to 
the body surface in the forebody regions was on the order of half of 
the local mean free path. The cell dimension normal to the surface 
in the wake regions was less than the local mean free path. Further, 
study of grid resolution is given in Appendix A of Ref. 19. Steady 
state was assumed when the number of simulated molecules in each 
region achieved a fixed value within fluctuations of 1 %. The average 
number of simulated molecules in a cell was 20 at steady state. The 
final results were obtained through a time-averaged solution over a 
large number of time steps. 

The molecular collisions are simulated by the variable hard sphere 
(VHS) molecular model. This model employs the simple hard- 
sphere angular-scattering law so that all directions are equally pos- 
sible for the postcollision velocity in the center-of-mass frame of 
reference. However, the collision cross section is a function of the 
relative energy in the collision. The freestream viscosity and mean 
free path are evaluated based on N 2 species using the VHS col- 
lision model with T re{ = 2880 K, d xt f = 3.08 x 10 -10 m, and 
s — 0.73. Energy exchange between the translational and inter- 
nal modes is modeled by the Larsen-Borgnakke statistical model 20 
with rotational and vibrational relaxation numbers of 5 and 50, 
respectively. 


Continuum Method 

The Navier-Stokes solver consists of an axisymmetric three- 
temperature, five-species implicit code. The set of equations solved 
consist of global mass, species mass, axial and radial momentum, 
global energy, rotational, and vibrational energy conservation equa- 
tions. The code uses Roe’s upwind scheme for the spatial inviscid 
fluxes, achieving second-order accuracy with the MUSCL scheme. 
The time integration is accomplished using an implicit lower- 
upper symmetric Gauss-Seidel (LU-SGS) scheme. This scheme was 
chosen because it only requires the inversion of a diagonal matrix 
at each point in the flowfield. With the large number of equations 
being solved, this is an attractive feature. A complete description of 
this method and modeling is presented in Ref. 15. 

To account for the effects of the Knudsen layer that develops 
on the surfaces in low-density flows, slip boundary conditions are 
used. The modeling employed accounts for velocity, temperature, 
and species concentration slip effects. A more detailed description 
of the slip boundary conditions used is given in Ref. 16. 

As shown in some previous comparisons between solutions using 
the DSMC method and the Navier-Stokes equations, 16 a meaning- 
ful comparison can be achieved provided the physical modelings 
employed by the two methods are matched as closely as possible. 
By doing this, the effects of using different physical models can be 
deduced, and the differences between the two solution procedures 
can be better observed. Thus, the chemical kinetic model and relax- 
ation rates employed in the DSMC method 17 are also employed in 
the continuum calculations. 


Freestream and Wall Conditions 

The freestream conditions considered are those experienced by a 
typical space vehicle at an altitude of 85 km during Earth entry and 
are listed in Table 1. Three configurations — 1) a 70-deg blunt cone, 
2) a 70-deg blunt cone with cylindrical afterbody, and 3) a proposed 
Mars Environmental Survey (MESUR) Pathfinder probe body — are 
considered in the present calculations. The base diameter is 2.0 m for 
each of these configurations. Therefore, the forebody configuration 
is the same as that proposed 21 for the MESUR Network probes 
(16 total). Since the current MESUR Pathfinder probe has a 2.65- 
m base diameter, the third configuration is the scaled Pathfinder 
probe. These configurations are shown in Fig. 1. The freestream 
parameters along with selected results are summarized in Table 2. 
The freestream mean free path is based on the VHS model, and it is 
calculated from the relation (see Ref. 13) 


(Too/r ref r 

[\/2n 00 o' re f(2 - co)°>r(2 - <z>)] 

Q) = S — \ 


( 1 ) 

( 2 ) 



R b =1 

.0 m 

R n /R b 

= 0.5 

R c/ R b 

= 0.05 

U R b = 

= 0.417 

R p/ R b 

= 0.5 

L p / R b 

= 1.0 

R ir/ R b 

= 0.215 

L m /R b 

m b 

= 1.057 


Fig. 1 Blunt cone and afterbody configurations. 


Table 1 Freestream conditions 3 


Altitude, 

Poo x 10 6 , 

Foo, 

TqOj 




Z-oo> 

km 

kg/m 3 

km/s 

K 

Moo 

Xo 2 

*n 2 

mm 

85 

7.955 

7.0 

180.65 

25.98 

0.2372 

0.7628 

7.189 


3 T w = 1000 K and Al = 28.96 g/mole. 
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a) Surface pressure 



b) Surface heating rate 



s/R n 

c) Skin friction 


Fig. 6 Surface results (altitude = 85 km and Ktioo = 0.00359). 


Comparison of DSMC and Navier-Stokes Results 

Comparisons of the chemical compositions from DSMC and 
Navier-Stokes calculations are shown in Fig. 2. Figures 2c and 
2d show that the Navier-Stokes calculations predict slightly lower 
values of atomic oxygen and atomic nitrogen in the wake re- 
gion compared with the DSMC results. The mole fraction of 
atomic oxygen in the wake region predicted by Navier-Stokes 
(with slip boundary condition) calculations is about 12% lower 
than that calculated by DSMC simulations (Fig. 7a), whereas 
the mole fraction of atomic nitrogen predicted by the Navier- 
Stokes calculations is about 18% lower near the base, and this 
difference decreases gradually in the downstream direction (Fig. 7b). 
The radial nondimensional density and axial velocity profiles 
(Figs. 4a and 4b) show very good agreement between DSMC 
and the Navier-Stokes calculations in the wake region. Figure 8 
shows that both the DSMC and the Navier-Stokes calculations 



a) Atomic oxygen profiles 



b) Atomic nitrogen profiles 

Fig. 7 Mole fraction profiles (altitude = 85 km and Ktioo ~ 0.00359). 



Fig. 8 Comparison of streamlines (altitude = 85 km and Ktioo = 
0.00359). 


predict a steady vortex in the wake region, and furthermore both 
calculations predict the same size (axial length), A = 0.77 (Fig. 5). 

Comparisons of DSMC and Navier-Stokes (with slip boundary 
condition) calculations of the surface pressure and surface heat- 
ing rate are shown in Figs. 6a and 6b. These surface quantities are 
shown in semilog scale as a function of nondimensional distance 
(s/R„) along the surface measured from the stagnation point to give 
better perspective of the magnitudes along the base plane. Along 
the forebody, Navier-Stokes and DSMC results are in good agree- 
ment for both pressure and heating rate. The surface pressure and 
heating rate along the base surface calculated by Navier-Stokes 
calculations are generally higher than the DSMC results (Figs. 6a 
and 6b). 



468 


DOGRA ET AL. : CHEMISTRY OF BLUNT-BODY WAKE 



Fig. 9 DSMC streamlines (altitude 



a) Pressure coefficient 



b) Heat transfer coefficient 



S/R n 

c) Skin-friction coefficient 

Fig. 10 Effect of afterbody configurations on surface results (altitude = 
85 km and Kn <*, = 0.00359). 



85 km and Kn^ = 0.00359). 


DSMC Calculations of the Blunt Body with Afterbody 
Configurations 

DSMC calculations were also performed for a 70-deg blunt cone 
with two afterbody configurations: 1) a cylindrical afterbody and 
2) a proposed afterbody configuration of the MESUR Pathfinder 
probe. These afterbody configurations were included to examine 
the impact on the surface heating and the flowfield features. 

Figures 9a and 9b show streamline plots for the cylindrical and 
MESUR afterbody configurations, respectively. It can be seen from 
these figures that the cylindrical afterbody has developed two steady 
vortices (Fig. 9a), whereas for the MESUR afterbody only one 
steady vortex exists in the wake as shown in Fig. 9b. The com- 
parisons of surface pressure, heat transfer, and skin-friction coef- 
ficients for the blunt body and two afterbody configurations are 
shown in Figs. lOa-lOc, respectively. The distributions of surface 
pressure and heat transfer coefficients on the MESUR afterbody 
surface are very close to the distributions on the base surface of the 
blunt body (Figs. 10a and 10b), whereas for the cylindrical after- 
body the heat transfer rate distribution on much of the afterbody 
surface is small compared with blunt-body base and MESUR after- 
body surfaces. 


Concluding Remarks 

Axisymmetric DSMC and Navier-Stokes (with slip boundary 
conditions) calculations of flows about a 70-deg blunt cone have 
been performed for flow conditions experienced by a space vehicle 
at an altitude of 85 km and a velocity of 7 km/s during Earth entry. 
Particular emphasis is given to the effects of chemistry on the near- 
wake structure and the surface quantities. The surface quantities of 
the afterbody surface are also analyzed. 

Essentially, all of the molecular oxygen and about half of the 
molecular nitrogen are dissociated before reaching the body surface. 
There is species separation among heavy (N 2 ) and light (O and N) 
species as the flow expands into the wake region. The density in 
the wake for the nonreacting gas model is slightly lower than the 
reacting model. The wake stagnation point for the reacting gas lies 
downstream of the nonreacting gas wake stagnation point. The base 
heating rate for the nonreacting gas is 240% larger than that of the 
reacting gas model with a noncatalytic surface. The comparison of 
the wake results predicted by the DSMC and the Navier-Stokes 
calculations shows good agreement for the flowfield quantities. The 
surface pressure and heating rate along the base surface calculated by 
the Navier-Stokes calculations are generally higher than the DSMC 
results. The presence of the afterbody has no effect on the forebody 
flowfield or the surface quantities. 
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